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Abstract 

We study the asymptotic behavior as n — > oo of the sequence 



1=0 

where B^^ and B^^ are two independent fractional Brownian motions, K is a kernel 
function and the bandwidth parameter a satisfies certain hypotheses in terms of Hi 
and i?2- Its limiting distribution is a mixed normal law involving the local time of the 
fractional Brownian motion B^^ . We use the techniques of the Malliavin calculus with 
respect to the fractional Brownian motion. 
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1 Introduction 

The motivation of our work comes from the econometric theory. Consider a regression 
model of the form 

yi = f{xi)+Ui, i>0 
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where (tij)i>o is the "error" and (xj)j>o is the regressor. The purpose is to estimate the 
function / based on the observation of the random variables yi, i > 0. The conventional 
kernel estimate of f{x) is 

X _ YJl=oKh{x^-x)y^ 

where K isa. nonnegative real kernel function satisfying K'^{y)dy = 1 and yK{y)dy = 
and Kh{s) = ■^i^(^). The bandwidth parameter h = hn satisfies /i™ — )• as n — )• oo. The 
asymptotic behavior of the estimator / is usually related to the behavior of the sequence 

n 

Vn = '^Kh{xi - x)Ui. 
1=1 

The limit in distribution as n — )• oo of the sequence 5„ has been widely studied in the 
literature in various situations. We refer, among others, to p^] and [9] for the case where xt 
is a recurrent Markov chain, to |151 for the case where xt is a partial sum of a general linear 
process, and [16] for a more general situation. See also |13j or |14j . An important assumption 
in the main part of the above references is the fact that n, is a martingale difference sequence. 
In our work we will consider the following situation: we assume that the regressor Xi = B^^ 
is a fractional Brownian motion (fBm) with Hurst parameter Hi G (0, 1) and the error is 
Ui = — B^^ where B^'^ is a fBm with H2 £ (0, 1) and it is independent from B^^. 
In this case, our error process has no semimartingale property. We will also set /i„ = n~" 
with a > 0. A supplementary assumption on a will be imposed later in terms of the Hurst 
parameters Hi and H2. The sequence Vn can be now written as 

n 

Snix) = Kin'^iB^^ - x)) [Bf_^\ - B^^) . (1) 

j=0 

Our purpose is to give an approach based on stochastic calculus for this asymptotic theory. 
Recently, the stochastic integration with respect to the fractional Brownian motion has 
been widely studied. Various types of stochastic integrals, based on Malliavin calculus. 
Wick products or rough path theory have been introduced and change of variables formulas 
have been derived. We will use all these different techniques in our work. The general idea 
is as follows. Suppose that x = 0. We will first observe that the asymptotic behavior of the 
sequence Sn will be given by the sum 

n 

an=Y. K{n'^Bl'^)K{n^Bf^)B [{Bf^\ - Bf-){Bf^, - Bf^)) . (2) 

i,jr=0 

This is easy to understand since the conditional distribution of Sn given B^^ is given by 

(a„)5Z 

where Z is a standard normal random variable. The double sum an can be decomposed into 
two parts: a "diagonal" part given by X^ILi -^^('^"-^2^^) ^ "non-diagonal" part given 
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by the terms with i ^ j. We wih restrict ourselves to the situation where the diagonal part 
is dominant (in a sense that will be defined later) with respect to the non-diagonal part. 
This will imply a certain assumption on the bandwidth parameter a in terms of Hi and 
i?2- We will therefore need to study the asymptotic behavior of 

n 

{S)n:=Y.K\n^B^^). (3) 
1=1 

(In the case H2 = ^ this is actually the bracket of S„ which is a martingale; this motivates 
our choice of notation.) We will assume that the kernel K is the standard Gaussian kernel 




This choice is motivated by the fact that K'^(n"B^^) can be decomposed into an orthogonal 
sum of multiple Wiener- Ito integrals (see [11], [3], [1]) and the Malliavin calculus can be 
used to treat the convergence of ([3]). Its limit in distribution will be after normalization 
the local time of the fractional Brownian motion denoted cL^^(l,0), where c is positive 
constant. Consequently, we will find that the (renormalized) sequence 5„ converges in law 
to a mixed normal random variable cWj^hj where is a Brownian motion independent 
from B^^ and c is a positive constant. The result is in concordance with the papers [15], |16j . 

But we also prove a stronger result: we show that the vector (5„, {Gt)t>o) converges in 
the sense of finite dimensional distributions to the vector (cW^. ^1(1,0)' iGt)t>o)i where c is 
a positive constant, for any stochastic process {Gt)t>o independent from B^^ and adapted 
to the filtration generated by B^^ which satisfies some regularity properties in terms of the 
Malliavin calculus. We will say that Sn converges stably to its limit. To prove this stable 
convergence we will express Sn 3S db stochastic integral with respect to B^'^ and we will use 
the techniques of the Malliavin calculus. We will limit ourselves in this last section to the 
case H2> \- 

We also mention that, although the error process B^"^ does not appear in the limit of 
(dl, it governs the behavior of this sequence. Indeed, the parameter H2 is involved in the 
renormalization of ([T]) and the stochastic calculus with respect to B^'^ is crucial in the proof 
of our main results. 

We have organized our paper as follows: Section 2 contains the notations, definitions and 
results from the stochastic calculus that will be needed throughout our paper. In Section 
3 we will find the renormalization order of the sequence ([1]), while Section 4 contains the 
result on the convergence of the "bracket" ([3]). In Section 5 we will prove the limit theorem 
in distribution for S'„(0) and in Section 6 we will discuss the stable convergence of this 
sequence. 
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2 Preliminaries 



Here we describe the elements from stochastic analysis that we will need in the paper. 
Consider H a real separable Hilbert space and {B{ip), (p £ V.) an isonormal Gaussian process 
on a probability space ($7, A, P), that is a centered Gaussian family of random variables such 
that E (B (if) B (ip)) = (99, ■(/')■«• Denote by the multiple stochastic integral with respect 
to B (see [10]). This In is actually an isometry between the Hilbert space 7^®" (symmetric 
tensor product) equipped with the scaled norm -^|| • ||-^®n and the Wiener chaos of order 
n which is defined as the closed linear span of the random variables Hn{B[ip)) where 99 E 
^) ll'/^ll'H = 1 ^iid Hn is the Hermite polynomial of degree n > 1 

^n(x) = i^exp^-j — ^exp^--jj, ..M. 

The isometry of multiple integrals can be written as: for m, n positive integers, 

^ {In{f)Im{g)) = n\{f,g)y^^n H TU = H, 

B{lM)Im{9)) = if m^n. (4) 

It also holds that 

/n(/) = /n(/) 

where / denotes the symmetrization of / defined by /(xi, ...,Xx) = ^ I]ae5„ fi^^{i)^- ■ ■ > Xa{n))- 

We recall that any square integrable random variable which is measurable with respect 
to the (T-algebra generated by B can be expanded into an orthogonal sum of multiple 
stochastic integrals 

F=Y,Ufn) (5) 

n>0 

where /„ G 7^®" are (uniquely determined) symmetric functions and /o(/o) = E [F]. 
Let L be the Ornstein-Uhlenbeck operator 

LF = -Y,nln{fn) 

n>0 

if F is given by 

For p > 1 and a G M we introduce the Sobolev-Watanabe space D"'^ as the closure of 
the set of polynomial random variables with respect to the norm 

= ll(^ - -^)^IIlp(q) 
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where / represents the identity. We denote by D the Mahiavin derivative operator that acts 
on smooth functions of the form F = g{B((pi), . . . ,B{ipn)) (g is a smooth function with 
compact support and tpi £ Ti) 

i=l * 

The operator D is continuous from D°'^ into 0"^-*^'^ [Ti) . The adjoint of D is denoted by 
5 and is caUed the divergence (or Skorohod) integral. It is a continuous operator from 
jja.p ^-^^ -j^^Q |jja-i,p_ have the following duality relationship between D and 6 

B{F6{u)) = B{DF, u)n for every F smooth. (6) 

For adapted integrands, the divergence integral coincides with the classical Ito integral. We 
will use the notation 

r 

6{u) = / UgdEs- 
Jo 

Let li be a stochastic process having the chaotic decomposition Us = X^„>o Inifni', s)) where 
/n(') s) G 7^®" for every s. One can prove that u G Dom 6 if and only if G for 
every n > 0, and "^'^=0 In+i{fn) converges in L'^{Q). In this case, 

oo oo 

6{u)=J2ln+iifn) and E|,5(n)p = ^(n + l)! 

n=0 n=0 

In our work we will mainly consider divergence integrals with respect to a fractional Brow- 
nian motion. The fractional Brownian motion (5/^)te[o,T] with Hurst parameter H G (0, 1) 
is a centered Gaussian process starting from zero with covariance function 

:=i(t2H^^2//_|^_^|2H)^ s,tG [0,r]. 

In this case the space T-L-}{ is the canonical Hilbert space of the fractional Brownian motion 
which is defined as the closure of the linear space generated by the indicator functions 
{l[o,t]5^ £ [Oi^]} with respect to the scalar product 

{l[o,t],l[o,s])nH=R"it,s), s,te[0,T]. 



3 Renormalization of the sequence Sn 

As we mentioned in the introduction, we will assume throughout the paper that x = in 
©, then 

n-l 

Sn := S„(0) = K{n^Bf^){Bf^\ - i?f ^). (7) 

1=0 
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We compute in this part the L^-norm of Sn in order to renormalize it. We have 

' n-l 



\i,j=0 

/n-l 

= E[J2K\n'^Bf^){B^^\-B 



H2^ 

j ' 

i,j=0 

n-l \ 



\i=0 / 

(n-l 
^ Kin-B^^)Kin-Bf^){Bf^\ - - i?f ) 

The summand T' will be called the "diagonal" term while the summand T" will be called 
"the non-diagonal" term. We will analyze each of them separately. Concerning T' we have 

Lemma 1 As n ^ +00, 

Proof: Through the independence off^fM and(5,^M 

V /t>0 \ J t>Q 

n-l 



T' = ^ E (i^^K^f )) E [{B^^\ - B^^f) . 

1=0 

SinceE((5,^^,-Sf^)2)=l, 

n— 1 n— 1 / -. 



where Z is a standard normal random variable. Recall that, if Z is a standard normal 
random variable, and if 1 + 2c > 

consequently, 

n-l 



27r\/l + 2n2"z2J^i 
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As n — 7- +00, T' behaves as such 

1=0 1=0 j=0 • 

-a-i/i+l /"l „-«-Hi + l 



n ^ ' TT ^ n 



27r^/2 Jo 27r^/2(l - ifi) 

The sign "~" means that the left-hand side and the right-hand side have the same hmit as 
n — )• +00. We win use this notation throughout the paper. ■ 

We will now compute the term T" . To do so, we will need the following Lemma (lemma 
3.1 p. 122 in [17J). 

Lemma 2 For every s,r £ [0,T], s >r and < H < 1 we have 

where fi = E{BfB^) and t > is a constant. 

Concerning the non-diagonal term of E (5^) the following holds 

Lemma 3 Suppose that 

a-AH2 + Hi + 2>Q. (11) 

Then, as n ^ +00, 

' t>0 V /i>0 



Proof: Using again the independence of (^Bf^^ and (^B^-^ 

n-l 

T" = ^B[K{n^B^^)K{n^Bf^))B[{Bll\-Bf-)iBf^,-Bf-: 

^ n— 1 

= ^J2^{Kin^B^^)Kin^Bf^)) fH,ii,j) 



where 



MiJ) = l\\i-j + M'""' + M'""' - 2 |i - j^'] . (13) 
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We need to evaluate the expectation E (K{n'^B^'-)K{n'^Bf^)) . Let T 



i'^' R{i,f 

,„^/iV wp hfl.vfi in 



be the covariance matrix of ( B^^ , B^'^) . We have |r| = (ij)^^^ — B?{i,j) and F ^ 



iri \-R{i,j) 



^ ^ '^afe'^^)- The density of (^Sf\Sfi) is then 

27r\/|r| 
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We obtain 



(2vr)V|r| 



e 2 e 2 e 



dxdy 



1 


V-^'v V 




1 




(2vr)2v 




1 








1 




(2vr)2v 




1 





e 2— e ^irr 



-— |_n +-[rr 



^.2 „2« + 2ffiUMM) 



dxdy 



^2 

n2«|r|+j2»l ^ 



dxdy 



R''{i,i) ..2 



2^irjTF^V (n2«|r|+j2-"l)2^' 



dxdy 



(27r)Vir| 

1 







2 






(„2. + i^jVp 







,2^1 



„2Qf|r|+j2Hi y 



dxdy 



(27r) 5 v^n2"|r| + j2Hi 



|rl(n2<^|r|+i2Wi) 



|r|v/(n2"|r|+i2/^i; 



27rVn2" |r| + j2-f^i v/(n2" |r| + i2^^i)(n2" |r| + p^^) - R^{i,j) 



Thus 



27rv/(n2° |r| + 22Hi)(^2a |r| + j2H^) _ i?2(^^ j-) 
1 

27rv^n4" |r| + n'^aj2Hi _^ ^2a^2Hi _^ I ■ 

Suppose that i > j. We use Lemma [2] to bound |r| = ^2_ffi^-2//i _ ^2^^ ^-^ from below. 
Therefore 



E 



(^K{n''Bf^)K{n'^Bf')^ < 



Since a2 + ^2 > 2a6 with a2 = n^°r(i - jf^^f^^ and ^2 = n'^'^{i^^^ + j2^^i) 

1 



E 



27rA/2V^n2"(i - j^i ^n'^'^if^^ + j2^^i; 
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and using the same inequality as above for = i'^^^ and 6^ = j"^^^ 



3a 



E(K{n"Bf^)K{n»Bf^)) < ^ \ n^ - (15) 

27rv2r4 (i — j) 2 j 4 i 4 

Since fH2{hj) behaves as H2{2H2 — — when z — j — )• c«, we can assert that 

2 V / 

Using p^ . we can write 

n— 1 n— 1 _3a 

2 ,. ..2H2-2 



i>j 27r\/2T4(i - j) 2 j 4 i 4 



and consequently 



2/^2-^-2 



n— 1 I ^ 

^// ^ H2(2H2 - 1) 3a r,rr ^1 o 3ffi , 1 \ n ) /A * 



— > C{Hi,H2)>0 



< 



H2{2H2-l)C{Hi,H2) _3^+2H2-^ 



47r\/2T4 

It follows that under condition (jlip ^x'^+^i-ij"" converges to zero as n — )• oo. 

As a consequence of Lemmas [T] and [3] we obtain the following L^- norm estimate for S'„. 
Proposition 1 Suppose that condition [Tl\) holds. Then, as n ^ oo 

1 



n 



27r^/2(l - Fi) ■ 

The condition (jlip will be discussed more thoroughly later (Remark 1, Section 5). 

4 The limit in distribution of {S)n 

Proposition [T] implies that the diagonal part of 5,^ is dominant in relation to the non- 
diagonal part, in the sense that this diagonal part is responsable for the renormalization 
order of which is n°'~^^^~^. As a consequence we need to study the limit distribution 
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of n"+^i-i(S')„ = n°'+^^-^^'^lQ K^in^Bf^). Using the self-similarity property of the 
fractional Brownian motion we have 



n— 1 rt— 1 

n 

i=0 i=0 



The limit of the above sequence is hnked to the local time of the fractional Brownian motion 
B^^ . For any t > and x G M we define L^^ {t, x) as the density of the occupation measure 
(see [1], m 

Jo 

The local time L^^ {t, x) satisfies the occupation time formula 

/ fiB^')ds= I L'''{t,x)f{x)dx (17) 
Jo Jm. 

for any measurable function /. The local time is Holder continuous with respect to t and 
with respect to x (for the sake of completeness L^^ {t, x) has Holder continuous paths of 
order 6 < 1 — H in time and of order 7 < in the space variable (see Table 2 in [5])). 
Moreover, it admits a bicontinuous version with respect to {t,x). 

Below, we give an important convergence result that will be necessary in proving the main 
result of this section. 

Proposition 2 The following convergence in distribution result holds 

Proof: Fix e > 0. Let Pe{x) be the Gaussian kernel with variance e > defined by 
Pe{x) = fY^ ^ • ^0^6 that for every s > 



V^n"+^^K2(n"+^^i?f^) = \v . (i?f^). (19) 
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Using (flUj) . we can write the left-hand side of as 



^ n-1 

j = n 

n-1 i+1 , X 

Tl — 1 i + 1 / \ 



j=0 

We win now estimate the three terms above and we wih show that each of them converges 
to zero (in some sense). Let us first handle the term We have 

Jo 2 JO 
It follows from [llj or [4J that 

[' MB^')ds ^e^o C 5o{Bf')ds = L''^{l,Q) (20) 
Jo Jo 

in L'^{Q) and almost surely, where L^^{1,0) is the local time of the fractional Brownian 

motion. Therefc 

be expressed as 



motion. Therefore an'^ clearly converges to zero as e — >• and n — )> 00. The term can 



- ) - j\e{B?')d^ (21) 

and for every e > it converges almost surely to zero as n — )> 00 using the Riemann sum 

(3) 

convergence. Let us now handle the term a\ given by 



= -Y.{p^^B'i")-nn~-i-^-ABl')] . (22) 

We will treat this term by using the chaos decomposition of the Gaussian kernel applied 
to random variables in the first Wiener chaos. Recall that (see [3], [6], [7], [12j ) for every 
G T~iHi (^Hi is the canonical Hilbert space associated with the Gaussian process 

PeiB^'K^)) = ^ Cmhm ^ -J (23) 

^ — ^ / ^ m-\-- 

m>0 I llifW^^ +e 



2 
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where C„ 



27r2'"m! ' 



Using this chaos decomposition, we can write Pe{B^^) — pi„-2{a+Hi) (S^^) as 



m>0 



[0,- 



Ey^ T / i(g)2m 



m>0 



[0,- 



n 



d'i,e,n' 



where 



H,e,n,m 



n) ' 



We win show that converges to zero in L'^{Q) as n — )• oo and e — )• 0. From (j22p one 

(3) 

can easily see that the diagonal part of converges to zero. We can also see, from the 

(3) 

expression of a„ , that the summands with j = vanish. Then, by using the orthogonality 
of multiple stochastic integrals ([10]), we obtain 

n-l 



m>0 

We can also write 



-2/fi(m+i) / -N -2Hi(m+i) 



di^e,n,mdj^e,n,m' 



E(a(3))2 



where 



n-l 



n 



m>0 



^Cl{2my.A^{e,n). 



m>0 



1 



n-l 



Am{e,n) = V ( -, - 



2m / . X -2/fi(m+|) / ■ x -2Hi{m+\) 



We can now claim that, for every fixed m > 1 



di^e,n,mdj^£^n,'i 



lim lim Am{£,n) = 0. 

e— >0 71— >oo 



(24) 



12 



Indeed, for every m > 0, we get 



V 



< 



< 



,2Hi 



.2Hi 



+ e 



m+1 



1 + 1 



+ 



(^) 



2Hi 



+ ln-2{a+Hi) 



2Hi 



+ 



,2H, 



,2Hi 



+ e 



+ 



n 



-2{a+Hi) 



Now, for every i,n,m, we have lim 

£-S>0 



and for every i > 1, 



n 



-2{cx+Hi) 



because a > 0. 



< 



n 



^2{a+Hi) 



m+1 



< C- 



n 



n 



2{a+2Hi) 



Furthermore, we know that 

n-l 



1 



i,i=0 



X 2m , 



n n 



n 



-2Hi(m+'^) ^ .~^-2H^{m+\) 



n 



converges as n — >• oo to /J /(J i?(u,f)2™(ui;)-2Hi(m+i)^^^^_ g-j^^g ^jj-g quantity is finite ([3] 
and [4]), it imphes (f24|l . 



We will now prove that 



^ C^(2m)!sup |^m(e,n)| < oo. 



(25) 



m>0 
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Relation ()24p and (j25p will imply the convergence of to zero in L^(r2). We need to find 
an upper bound for the terms \di^s,n,m\ and \dj^£^n,m\ in order to continue. 



,2Hi \ 2 / I i\'iHi 



One can note that 



(1 + en?^i-'^H) J I (1 + in-2«i-2^^) 



< ( ]„. pffj ' < 1 and < ( ) < 1 

because en?^i~'^^ > 0. From the above inequalities, we can deduce that 

- V (1 + en2^^i-2^^) ) (1 + l„-2a^-2i/) j - 

and finally, 



|'^i,e,n,m| ^ 1 and |dj,e,n,m| ^ !• 



By bounding from above the terms and by 1 in X]„>o C^(2m)! sup„ 

we obtain that 

/• • \ 2m /•\-2Hi(m+i) / • \ -2_f/i(m+i) 



5^C^(2m)!sup|A^(e,n)| < Ci(2m)! sup E ^^4^'^ U U 

m>0 "'^ m>0 " j,i>l,j7^i V / V / V / 

n~l , / ' \ \ 2"i / ■ \ —2H\m / ■ ■ \ ■ 

m>0 *J>l,j7^i VV// V/ V/ 

Let's focus on the case where ifi < ^ first. Let Qhi {z) be the function defined by 



Qhi {z) 

For i^i < i, we have 



if ze(o,i] 

if z = 0. 



Qh, iz) < z'^K 
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Indeed, the function f{z) = 1 — z^^^ — (1 — z)^^'^ is negative on [0, 1], increasing on [i, l] 
decreasing on [O, ^] and /(I) = /(O) = 0. It follows that 



-, n-l ^ . X 2Him / ■ • \ -Hi 

Y,Cii2my.sup\Am{e,n)\ < 2 J] C^(2m)! sup - J] K 

-"^n ' m>0 i,j=0]i>j ^ ^ ^ ^ 

1 /■\Hi{2m-l) /■\- 

m>0 i,j=0;i>j ^ ^ ^ ^ 

m>0 i=0 j=l 

< 2 VC^(2m)!supn2^i-2 Vr^i(2-+i) f xi^^i^^-^)dx 

< 2^C^(2m)!sup- 

m>0 

< 2^C^(2m)!sup. 



m>0 ' m>0 i,j=0,i>j 

1 /„-\^^l(2m-l) -Hi(2m+1) 



m>0 i=0 



m>0 - _L _L . 

n-l 



2Him - 




+ 1 


n~ 


1 




2Him- 


Hi 


+ 1 


1 






2Him- 


Hi 


+ 1 



m>0 ± ± ■ -^Q 



Cl{2m)\ 



,„ . _ ^2Him-Hi + l 

m>0 ^ ^ j^>Q 



Given that, by using Stirling's formula, the coefficient C^(2m)! behaves as we obtain 
that the above sum is finite. Thus, we obtain the convergence of a^f^ to zero in L^(r2) for 



Hi < 



Let us now treat the case ffi > |. We know (see [4], Lemma 1) that the function Qh 
is increasing on [0, 1]. Since 2 < = 1 - i it holds that Qh(^) < Qh{1 - i)- Then 

^ n-l / n-l / . . N -Hi 

Y,Cl{2m)\^MAm{e,n)\ < 2 J] C^(2m)! sup ^ 1 - - E 

m>0 m>0 1=1 j = l 

= 2j:ci(2™)!.upij:<3„(i-i)(i) E L 

m>0 i=l ^ / \ / j=i-' — 

< CHY^Cl{2m)\snp-Y,QHll--)"^ 

m>0 t=l 

1 "^^ / 1\ / i ^ 1-2//1 
~ CHVC^(2m)!sup-J^gH 1-- 

\r, n \ I J 



m>0 j=l 
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By adapting Lemma 2 in [J] (by separating the sum over i in a sum with i < (5 and ^ > S 
with 5 suitably chosen), we can prove that 

-I n-l / / -x \ 2m / -x —2Him / ■ -x -H 



n ^ — ' \ \t J J \t J \nn 

i,j=o \ \ / / \ / \ . 

with c{Hi) not depending on m nor n. As a consequence 

V c^(2m)!sup|A^(e,n)| < c(Fi)c^(2m)!m"^ . 

m>0 

The Stirhng formula implies again that the above series is finite. ■ 

Theorem 1 Let {S)n be given by Then, as n ^ oo, we have the convergence in 

distribution 

n^+^^-\S)n^ [ K\y)dyL^^ {1,0) 

JR 

where L^^{1,0) is the local time of the fractional Brownian motion 

Proof: Using Proposition [2] it suffices to check that n°'+^^ K'^{n°'~^^^B^^)ds converges 
to J^K'^{y)dyL^^ (1,0). Using the occupation time formula (fT7|l . we obtain 

f\^{n'^+HiB^i)ds = n'^^^^ [ K^{n''+^^x)L^^{l,x)dx = [ K^{y)L{l,yn-''-"^)dy 
Jo Jm. Jr 

which converges as n — )■ oo to K'^{y)dyL^^ (1, 0) by using the continuity properties of the 
local time. ■ 

5 Limit distribution of Sn 

In this paragraph, we prove the limit in distribution of ([7]). Recall the notation (113p and 
let's consider the Gaussian vector 

= = {B^'-B^\...,B^^-B^l,). 

From this definition, it follows that 

n—l n—l 

Sn = Y. K{n^B^^){Bf^\ - Bf-) = ^ K{n^B^^)X^_^\. 

i=0 1=0 
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Theorem 2 Let (Sn) be given by ^ and assume that 

a<l-Hi (26) 

Then we have the convergence in law 



n 



n—^+oo 



(1,0) 



where L^^{1,0) is the local time of B^^, di := J^K'^{y)dy and W is a Brownian motion 
independent from B^^ . 

Proof: We will study the characteristic function of nz"*" 2 2 5'„. In order to simphfy the 
presentation, we wih use the following notation. Let be the imaginary unit and A„ be 

A„ = with A G M. 

Using the independence of the two fBms and computing the conditional expectation of 



because if X is a Gaussian vector with mean /i and covariance matrix S, it's characteristic 
function is given by 



It follows that, with /i^jl^jj) given by p^ . 



Consider the process {Vn)n>o defined by 

K (""^w ) ^ ) 1^ - ^1'^'"' ^^^^ 

and the function ■i/^ defined by ■0(3^) = e"'*'"''^^*^^''^^""'^^^. Note that, since we excluded the 
diagonal, the integral duds in the expression of y„ makes sense even for H2 < ^. Note also 
that is a bounded variation process (its quadratic variation is 0). Furthermore, 
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Using the change of variables formula for bounded variation processes, it follows that 

V'(K) = 1 + / i^'{Vs)dVs 

Jo 

i.e., 

-XlH2{2H2-l)V„ 



e 



Jo 



Therefore 
E 



= E ( e 2 



E(Ti)-E(r2). 



We will now focus on the term E {T2) and show that 



rl 



From 



dVs = K (n"5[^^) K [n'^Bfj^) \s - j ds 



we get 



E(r2) = ^(^\lH^(2H2-l)e~^~^^^=°^'(''^^^') 

X y e-^"^^(2i/,-i)y. y ^ ) 1^ - ^r^"' ^^^^^^ j 

Recall that the following holds 
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This can be seen for s integer as at the beginning of this proof and also (|27p can easily be 
checked for any s > 0. We will use this property to compute the following upper bound for 
E(|T2|) 



EdTal) < E 



<1 



X I k( n'^Bfn K [n'^B^A H2\2H2 - l\\s - u\^^^-^ duds 



n / \ 



X i'K( n'^B^J ] K (n'^BfA H2\2H2-l\\s- u\^^^~^ duds 



This is true because all the terms of the double integral are measurable with respect to the 
filtration generated by {B^^,u < s). At this point, we use to write 



\s — ul"^^^ ^ duds 



< E 



< E A 



K ( n^Bf^' ] K ( n'^Bf} ] Ho \2H2 - Ills - ul^"^'^ duds 



2H2-2 



n [■\s\ 



10 JO 

n— 1 i— 1 



K ( n°-B^i ) K ( n'^B"^ ] Ho \2Ho - 1| |s - u 



\2H2-2 



duds 



< BixlY,Y.K(n''B^'^K(n''Bf'^H2\2H2-l\J^ J \s - ul"^^''^ duds 



1=0 j=0 



Assume that H2 > ^, ergo \2H2 — 1| > and fH2{hj) > 0. Consequently, 
E(|T2|) < E j^j;j;K(n-i3f^)i^(n"Sf 

n— 1 i— 1 



i=0 j=0 



< E I yn"+^i- 



i=0 j=0 
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The previous term is exactly the non-diagonal term of the L -norm of n 2 ^ 2 2 5„ and we 



know that under condition (|lip . it converges to zero when n — )• +00. Finally we have 

E(|T2|) 0. 

Assume now that H2 < \. It follows that \2H2 — 1| < and fH2{hj) < 0; which gives us 



i=0 j=0 



■2 n-1 i-1 



\ i=o i=o 

As in the previous case, this term is again exactly the non-diagonal term of the L^-norm of 

a I -f^l 1 

n^^~~2 Sn and for the same reasons, we get the following result again (which is now valid 
for any H2 G (0,1)) 

EdTal) —> 0. 

n— >+oo 

Concerning the term Ti, we note that 

E(ri) = E ^e^^^-^^' 
and the result follows from Theorem [TJ 



Remark 1 The following comments deal with the conditions Ul\} and i26\). Condition 
i26\) is a natural extension of the condition a < | in e.g. J 15^ - U6f which means that the 
bandwidth parameter satisfies nh'^ = nn~'^" — ?• 00 as n —)■ 00. From and i26\) . this is 
the constraint we find for a (considering a is our degree of freedom) 

0<Hi<l 
0<H2<1 
a>AH2-Hi-2 
a <l- Hi 

As an example, consider the case where H\ = H2 = H. Those constraints become 

r o<H <i 

\ {3H -2)^ <a<l-H. 
For this system to have a solution, we need to verify that 

3H -2 <l- H ^ H <-. 

4 



< i^i < 1 
0<H2<1 
AH2- Hi-2 < a <l- Hi. 
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As a result, our constraints become 



< < I 
(3i/ - 2)+ < a < 1 - H. 

We could also consider the case where a has a fixed value and where the constraints would 
be on Hi and H2 ■ 



6 The stable convergence 

In this section we will study the convergence of the vector {Sn, {Gt)t>o) where {Gt)t>o is a 
stochastic process independent from B^^ and satisfies some additional conditions. In this 
case, since the process {Gt)t>o is not necessarily a Gaussian process and since no information 
is available on the correlation between B^^ and Gt, the characteristic function of the vector 
{Sn, {Gt)t>o) cannot be computed directly. To compute it, we will use the tools of the 
stochastic calculus with respect to the fractional Brownian motion. The basic observation 
is that Sn can be expressed as a stochastic integral with respect to B^'^. Indeed, 

n— 1 n— 1 

Sn = ^KKSfO«\-i?f^)=^i^KSf^)<^^^(l[,,+l](.)) (28) 
1=0 i=0 

= X:<^''^(^K5f^)l[^.+i](-)) + ( zj^^i^(n"i?f^) ,lt,,+i](-)^ 



0fromB»i±Bf2 



n—1 n—1 



^ / K{n^Bf^)dB^-=Y^ / Kin^Bf^B^- = / i^(n«Rf i)di?f ^ 



.-O"' i=0-"' 

We will also use the "bracket" of Sn- This quantity equals 

n-1 n-1 ^j+i 



Y^KHn-Bfn = . K\n-Bf^)ds 

i=o 1=0 

ri+l rn 

= Ey. K'in^Bf^Ods = K\n'^Bf^l)ds. 



Before going any further, we will describe the elements of the stochastic calculus with 
respect to fractional Brownian motion that wc will be using in the sequel. We will start by 
introducing some notations and definitions. Let (p be the function defined by 

(t){s,t) = H{2H -l)\s-tf"~^ , s,teR. 
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Let D (introduced in section 2) be the Malliavin derivative operator with respect to the 
fractional Brownian motion with Hurst parameter H. Based on this operator, let D'f' be 
another derivative operator (called the (/)-derivative operator) defined by 



(j){t,v)Dj,Fdv 



for any F in the domain of D. For more details about this operator, see [2]. Let £^(0, T) 
be the family of stochastic processes F on [0, T] with the following properties: F € /^^(O, T) 



if and only if E 



\F\ 



and E 



rT rT 
Jo JO 



DfF, 



< CO, F is (/>-differentiable, the trace of DtFt,0 < s,t < T, exists. 



dsdt 



< oo and for each sequence of partitions (tt^, n G N) such that 



as n 



+00, 
n-l 

EE 

i,j=0 



IX IX DtFJ^^DtFJ^^-DtFtDtF, 



dsdt 



and 



tend to as n 



E 



IF" 



+00, where 7r„ :< = 4"^ < 4"^ < - < 4-4 < 4"^ = T. 



In our particular situation, we are dealing with processes of the form FudB^ + Gudu, 
(where is a fractional Brownian motion with Hurst parameter H), for which the follow- 
ing Ito formula holds in the case H > ^. 



Theorem 3 Let r]t = /g FudB^ + Gudu, for t e [0, T] with E 



sup I Gs 

0<s<T 



< OO and let 



{Fu,0 <u<T) he a stochastic process in C^{0,T). Assume that there is a f3 > 1 ~ H such 
that 



E 



\F — F 



<G\u-v 



2/3 



(29) 



0. 



(30) 



where \u — v\ < (" for some C > and 

Urn E Dt{Fu-F,] 

0<u,v<t,\u-v\^0 I 

Let f : M+ x M — t- M be a function having the first continous derivative in its first variable 
and the second continous derivative in its second variable. Assume that these derivatives 
are bounded. Moreover, it is assumed that E 

[0,r]) g£<^(0,T). Then forte [0,T], 

ds ' / Qx 



fit,Vt) 



/(0,0) + 
rt Qf 



f 

Jo 

-{s,ris)ds + 



FsDtv. 
'df 



ds 





+ i + I g^is,Vs)FsDUds. 



< 00 and {^^Fs,s G 
{s,ris)Gsds 



(31) 
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We also have the fohowing technical lemma P-71) that will be particularly useful for 
our future computations. 

Lemma 4 Let {Ft,t G [0, T]) be a stochastic process in £,^(0, T) and 

2" 



sup E 



0<s<T 



< DO 



an 



d let r]t = /o FudBll for t G [0, T\ . Then, for s, t e [0, T] , 

Dfrjt = f DfFudB^ + /* Fu<P{s, u)du. (32) 
Jq Jo 

It is now possible to state the main result of this section 

Theorem 4 Assume that and i26\} holds. Let {Gt)t>o be a stochastic process inde- 
pendent from B^^ and adapted to the filtration generated by B^^ such that for every t > 
the random variable Gt belongs to B^'^ and ||Z)<jGi| < C for any s,t and to. Then the 
vector {Sn, iGt)t>o) converges in the sense of finite dimensional distributions to the vector 
{cW^Hii^i 0)' (G't)t>0; where c is a positive constant. 

Proof: In order to simplify the presentation, the following notations will be used. We will 
denote by A„ (like we did in a previous proof) the quantity 

a , Hi 1 

A„ = An 2 2 2 

where A G M. The following notation will also be used: 

e(A,n) = e-^""'''''"'^o"^'(""<] )'^" = e-4"'^+"^~' ^'(-"^f M. 

Let {Ft,t > 0) and {Gt,t > 0) be two stochastic processes defined by 

r F^ = Kin-Bf^l) 

\Gu = -io\K\n^Bfj^) 

and let (?7p"\t > 0) be the stochastic process defined by 



FudB^ + 1^ Gudu = K{n"Bf^^)dB^^ -io^ 1^ K\n''Bfj^)du. 

Consider the function / : C — ?■ C, /(x) = e*"'^"^. We can apply the Ito formula to firfj^"^^) 
in order to obtain 

\2 ft 



■ \ (An) 



1 + ^^ e^o^-^'"-' K\n^Bf^l)ds 
+ioXn e^«^"''^^"'i^(n°Sgi)dSf 2 

-Xl e^o^"''^^"^K(n"5gi)Z)f (33) 
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where D'f''^'^ is the operator D"^ introduced above with respect to the fractional Brownian 
motion 

_ We use Lemma m to compute Df'^'^'r]i'^"\ We get 



=0 from BHi±b"^ 



\2H2-2 



=0 from B^i±Bf2 



= H2{2H2-l) f^K{n^Bfj^)\s 
By substituting in (j33|) . we obtain 

e.oA„,^' = 1 + ^ [\^oxW-'KHn»Bfms 

-XlH2{2H2 - 1) ^* e^o^"^^'"'K(n°Bgi) 

X / i^(n"5jp|s-u|2^^~^dM(is. (34) 



By multiplying both sides of (j34p by e ^ = e(A,n), we obtain 

gioA„ /o" (n-B^f 2 ^ 4 /„" )du 

2 Jo 

f " (A ) 

+ioXn j e*"^""^ " K{n^B^^I )dB^^ • e(A, n) 
-XlH2{2H2 -I) J e'°^"^^ K(n"5gi) 

X / K{n''Bll^)\s-uf"^-^duds-e{X,n). (35) 
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The sum of the last two terms in (j35p can be written in a more suitable way by using sums 
instead of integrals. Together, these two last terms give us 



ds ■ e(A, n] 



,a T3H1 







du 



+ / Kln'^Bf'] \s - u\'^"^~^ du 



ds ■ e(A, n) 



i— 1 /"i+l 

xF2(2i72 - 1)E^ (""-^fO / \s - u\^"^~'^ duds ■ e{X,n) 

(n— 1 fi+l 

i-/72(2i/2-l)^%-n|' 



1 2/^2 -2 



ds ■ e(A, 



(36) 



Let us now fix . . . , /3iv e M and ti, . . . , ^at > 0. We need to show that E [e^o>^"^nfJ-o{PiGf,+-+PNGt^) 
converges to E e 2 g«oipi(Jtii----i-P]vLJtivi ) . "We wih use the notation 



gN ■= e 
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By combining relations (|35p and (|36|) . we can write 



E 



^Zq An 5"^, 



aN 



E(e(A,n)5Ar) 



-E ( ioA„ 

n-l 



\ 1=0 ■'^ 

i 1 

xH2{2H2-l)^K {n'^Bf'^ j \s - u\'^"^-^ duds x e{X,n)gN 



n-l 



j=0 
fi+l 



du 



ds X e(A, n)g]\f 



X ^-H2{2H2-1) J' Is-uf"" 
:= E(e(A,n)g^) + ^i^ + ^2* + ^3^ 
Let us begin by proving that the term T2 converges to zero as n — ^ 00. Since 



(37) 



(38) 



for every s < n and since |e*°^| = 1 for every x real, can be bounded as follows 

n-l 



< E A2c(i/2)EE^M""^fO^K^. 



\s -ul'^^^-'^duds 



n — 1 i— 1 
i=0 j=0 



i+1 n+i 



Is - ul'^^^-'^duds 



and this goes to zero as in the proof showing that the non-diagonal term goes to zero under 
the renormalization n'^'^^^^^. Let us now handle the term T^. By using the independence 
of and 

we can write 

= E (^ioXn £ e^«^"''^'"'i^(n-i?gi)e(A, n)dBf^ • g^^ . 

The duality formula is used to obtain 

= E(ioA„(l[o,„]e*°^"''''"'K(n"i?ff^)e(A,n),D^^<7^)^^^ 



E 



A„(l[o,„]e*"^"''-''"'K(n-i?f{^)e(A,n),57v 



N 



k=l 
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^ Jo Jo 



r — ul"^^^ ^ (jirtpududr 



Recall that the following formula holds (see [lU] for further details) 

r-T l-T 

'0 Jo 

for any pair of functions in the Hilbert space T-Lh2 ■ This formula is used to write T* as 

where the fact that Gt is adapted to the filtration of B^'^ is used. It suffices to show that 
for every fixed t > 0, 

converges to zero as n — ?■ oo. Since the derivative of Gt is bounded and using (j38p we find 
that the above term is less than 



XnCiE ^ - vl^'^'-^dvdu^ 

n—l „ 

A„ci^E(K(n"i3f^)) / 



i=0 

n-l 



= XnCiCH,Y.B[K{n"Bf^)j + 1 - f + \i - f + \i + Ip^^ - i^^^) 

1=0 

where ci is the constant upper bound of the derivative of Gt and cj^j is a constant depending 
only on H2. Since for every fixed t > the function (-|i + 1 - i^-^^^ + K - tp^^ + \i + Ip^' - i^^^) 
behaves, modulo a constant, as i^^^"^ and since the order of the expectation of K{n°'B^^) 
is the same as that of n~°^i~^^ it is clear that Ti converges to zero as n — )• 00. 

Finally, we will show that T3 converges to zero. Note that the term T3 can be expressed as 
follows 

= XlY^^(^K\n-B^^) f^\^^>^-^^''' {^--H2{s-^r^^^ 
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At this point, we will again apply the Ito formula for e "''^ . It implies that 



i=0 



2H2-1 



ds ■ e(A, n)gN 



(n—l / -1 



is-i) 



2H2-1 



K{n''Bfj.)dB^^ds-e{\,n)gM 



xH2{2H2-l) / e 



Kin'^Bf^^) I K{n^B^^^)\u-v\''''-''dvduds■e{Kn)g^ 



7H1 



2H2-2 



= 5« + 6(2) + 6(3) + 6(4). 
The first summand 6(^) vanishes because the integral 



vanishes. The second summand 6(2) goes to zero as n — )• 00 using exactly the same argument 
as for the convergence of T^. Concerning the third summand, b^^\ using (j38p and the fact 
that \gj\[\ = 1, we get 



6(3) < E 



< E 



XtY^K' (n"Sf^) |i - H2{s - ^r^-'\ [ ^^n^Bfj^ 



)duds 



<i 



1 /■«+! / 



i=0 

n— 1 i— 1 



(S - 

<1 



ds 



n— 1 i— 1 



< E I A^ ^ ^ in-Bf^ ) in-Bf ) I + E I A^ i^^ n"i?f ^ 



i=0 j=0 



j=0 j=0 
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The second term goes to zero because E E5=o^^ (™"^»^')) behaves as n ^Ji-r-L 

and the first term goes to zero because the non-diagonal term is dominated by the diagonal 
term. Analogously to the convergence of , the last summand in (I39p converges to zero. 
This completes the proof. ■ 
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